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Abstract 

The effect of hydrodynamic interactions on the free-solution electrophoresis of polyelec- 
trolytes is investigated with coarse-grained molecular dynamics simulations. By comparing 
the results to simulations with switched-off hydrodynamic interactions, we demonstrate their 
importance in modelling the experimentally observed behaviour. In order to quantify the hy- 
drodynamic interactions between the polyelectrolyte and the solution, we present a novel way 
to estimate its effective charge. We obtain an effective friction that is different from the hy- 
drodynamic friction obtained from diffusion measurements. This effective friction is used 
to explain the constant electrophoretic mobility for longer chains. To further emphasize the 
importance of hydrodynamic interactions, we apply the model to end-labeled free-solution 
electrophoresis. 
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1 Introduction 

Nowadays, electrophoretic separation methods are widely applied to study polyelectrolytes 
(PEs) such asproteins, DNA and synthetic polymers [Q], Q, [3Q- While there exist several the- 
ories 0, HI, [(I Q] that have been successfully used to describe qualitatively the experimentally 
observed behaviour of various PEs, there are still many open problems to address. 
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Recent experiments on strongly charged flexible PEs, such as polystyrene sulfonate (PSS) 
and single- stranded DNA (ssDNA) of well defined length have shown a characteristic be- 
haviour for the short chain free- solution mobility /x fill, 10, 11]: after an initial increase of 
the mobility with increasing length, /j passes through a maximum, and then decreases towards 
a constant mobility for long chains. 

The increase for short chains and the long-chain limit can be explained within the theoretical 
approaches, but the origins of the maximum for intermediate chains remain indistinct. To some 
extent this can be attributed to the simplifying assumptions made in those models regarding 
the interplay of the various interactions: the Coulomb interaction between the charged PE 
monomers and its counterions, the external electric field likewise acting on the charges, and 
the hydrodynamic interactions with the solvent. 

To provide a fundamental understanding of the involved dynamics, we believe it is manda- 
tory to study the effects of these forces on a microscopic level, thereby taking into account 
full electrostatic as well as hydrodynamic interactions. Continuing the work of [12], we use 
coarse-grained molecular dynamics simulations to determine the transport coefficients and 
structural properties of strong polyelectrolytes in free-solution electrophoresis. We charac- 
terize the hydrodynamic interactions between solvent and solute on a microscopic level, and 
determine the relevant length scale for these interactions. 

The article is structured as follows: in section [2] we introduce the simulation model and 
specify all relevant parameters. After comparing the simulation results to experimental mea- 
surements of the PSS transport coefficients for varying length (Section 13.11) . we switch off 
hydrodynamic interactions to characterize the behaviour in absence of them. In section 13.41 
we introduce the concept for an effective friction of the polyelectrolyte counterion compound 
with the surrounding solvent. To quantify the effective friction, we propose a novel way to 
estimate the effective charge. Finally, we show that the effective friction, measured in this 
way, deviates from the hydrodynamic friction that can be obtained from diffusion measure- 
ments, and that this deviation is the reason for the observable length independent mobility of 
long polyelectrolyte chains. Finally, in section [3771 we demonstrate, how increasing the ef- 
fective friction by a drag label can restore size dependent mobility. We finish this paper with 
concluding remarks. 



2 Simulation model 



We use a charged bead-spring model within the Espresso package |]13ll . to study the behaviour 
of flexible linear PEs of different lengths. All parameters are given in reduced units with 
energy scale kT = 1.0 and relevant length scale <To that is used to match the model to a 
specific polyelectrolyte. For this paper, we chose er — 2.5 A, which is the distance between 
two charges along a fully sulfonated PSS backbone. The simulation time step is r = 0.01. 

The PE is comprised of N negatively charged monomers carrying a charge of — le , where 
e is the elementary charge. 

The monomers are connected by finitely extensible nonlinear elastic (FENE) bonds £/fene = 
|/c.R 2 ln ^1 — (-^) 2 j, with stiffness k = 30, and maximum extension R = 1.5, and with r 
being the distance between the interacting monomers B14I1 . Additionally, a truncated Lennard- 
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Jones or Weeks-Chandler- Anderson (WCA) potential U u = 4e 2 - (^) 6 j, with depth 

e = 0.25 and width a = 1, is used for excluded volume interactions jl5\ \. 

For charge neutrality, iV monovalent counterions of charge +leo are added that are subject 
to the same LJ potential. All particles have the same size. 

The simulations are carried out under periodic boundary conditions in a rectangular simula- 
tion box. The size L of the box is varied to achieve a constant monomer concentration cpe and 
counterion concentration cci independent of chain length. To match experimental conditions 
with a PSS monomer concentration of 1 g/1 or 5 mM, we vary L between 34 (N = 2) and 89 
(N = 32). 

Full electrostatic interactions are calculated with the particle-particle-particle mesh (P3M) 
algorithm lllql . The Bjerrum length l B = el/AnekT = 2.84 in simulation units corresponds 
to 7. 1 A (the value for water at room temperature). 

To account for hydrodynamic interactions (HI), we frictionally couple all particles to a 
lattice Boltzmann fluid as detailed in QJ/ZD - The modeled fluid has a kinematic viscosity v = 
3.0, a fluid density p = 0.864, and is discretized by a grid with spacing a = 1.0. The coupling 
parameter is r bare = 20.0. 

In order to further characterize the impact of HI on the system's dynamics, we compare to 
simulations with a simple Langevin thermostat that does not recover long-range hydrodynamic 
interactions between the monomers, but only offers local interactions with the solvent. We set 
the friction parameter T = 15.34 to match the single particle mobility of the Langevin system 
to the one with full HI. 



2.1 Determining transport coefficients 

Using this model, we determine the single chain diffusion coefficient D and the electrophoretic 
mobility ji of the model PE from simulations without an applied external field using the fol- 
lowing Green-Kubo relations: 

oo 

1 



D = 3 / (v c (r) ■ v c (0))dr (1) 



o 



3k~r^ qt J ( v "c(r)-vl(0))dT (2) 







Here, v c is the center of mass velocity of the PE, vl is the velocity of every charged particle 
in the system, i.e. monomers of the PE and associated counterions, and qi is the corresponding 
charge. The ensemble average ((. . .)) is taken over 10 4 statistically independent samples. The 
integrations use analytic fits for the slowly converging long-time tails. 

The advantage of using a Green-Kubo relation (2) for the electrophoretic mobility is that we 
can obtain both transport coefficients from the same simulation trajectories at zero field. Do- 
ing so we avoid conformational changes to the chain structures or the counterion distributions 
by an artificially high external electric field, which is sometimes used in other simulations 
to separate the directed electrophoretic motion from the Brownian motion within reasonable 
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computing time. This approach was successfully applied in simulations to determine the elec- 
trophoretic mobility of charged colloids lll8L ll9n. 



3 Results 
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Figure 1: The normalized diffusion coefficient D/D for PSS of different lengths N as ob- 
tained by electrophoresis NMR agrees with the simulation results with full hydrody- 
namic interactions (HI). 



In an earlier publication this year [12], we compared the simulation results to experimental 
data for short PSS obtained by two different experimental methods, namely, capillary elec- 
trophoresis lIsiQ and pulsed field gradient or electrophoretic NMR [llE 0, S &4f]. 

In Figure [Q we compare the diffusion coefficient D determined by simulations to the 
experimental results. Both results are in good agreement and exhibit a power law scaling 
D = D N m , with an scaling exponent m = 0.63 ± 0.01, which agrees with previous results 



'Please refer to Jl2ll for details on the experimental conditions. 
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Figure 2: Including hydrodynamic interactions (HI), the normalized electrophoretic mobility 
/i//!fd as a function of the number of repeat units N obtained in simulations re- 
produce a maximum for intermediate chains as well as the long-chain behaviour 
observed in experiments. 
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fill I25L |26[ ll^]. Here, the simulated data is normalized by D = 0.052, the monomer dif- 
fusion as obtained by a power law fit, and the experimental data by the monomer diffusion 
coefficient of D = 5.7 ■ 10~ 10 m 2 /s. 

Figure |2] shows the electrophoretic mobility /i normalized by /^fd> the constant (not length- 
dependend) value for long chains. We compare the simulation results to two diffent experi- 
mental data sets. The simulation results reproduce a maximum for intermediate chains as well 
as the long-chain behaviour (the constant mobility for long chains) observed in experiments, 
if the hydrodynamic interactions are properly accounted for. 

"7. O A heonro r»f hurlmrlwnamir intorar^tinnc 
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Figure 3: The normalized diffusion coefficient D/D with and w ithout hydrodynamic inter- 
actions (HI). With HI the experimentally observed sea ling is recovered; without HI 
the chain diffusion is Rousse-like. 

Figure|3]illustrates the difference in the scaling of the diffusion coefficient with and without 
hydrodynamic interactions. For the latter, Rousse behaviour is observed with a scaling expo- 
nent of m = 1.02 ± 0.02 for the chain diffusion. Likewise, as shown in figure|H the neglect of 
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Figure 4: The normalized electrophoretic mobility without hydrodynamic interactions (HI) 
strongly deviates from the one obtained with HI. The plotted curve, equation ©, 
describes this observation by using the local force picture. 
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hydrodynamic interactions leads to a decreasing electrophoretic mobility for short chains and 
a constant value for long chains. The same observations were made in a recent publication by 
Frank and Winkler Q. 



3.3 Local force balance 

The behaviour in absence of hydrodynamic interactions can be understood in a local force 
balance picture. A strongly-charged polyelectrolyte such as PSS is surrounded by a cloud of 
oppositely charged counterions, some of which move with the chain, thus forming a molecular 
compound with the PE that has a reduced electric charge. This effect is also know as counte- 



rion condensation 112911 . Without hydrodynamic interactions, every particle of this compound 
is subject to the same frictional force specified by the Langevin thermostat. For a not to strong 
applied electrical field E, the electrophoretic mobility can be obtained by /i = v/E, where v 
is the steady-state velocity of the compound (PE and co-moving counterions). The two acting 
forces, the solvent drag force -F Solvent and the electric driving force F Fie i d , are of equal magni- 
tude and opposite direction. We now define iV CI to be the number of co-moving counterions. 
This leads to the following equation for the forces, where ±Q is the charge of the monomers 
and the counterions, respectively. 

^Solvent = -l>(iV + N C l) (3) 
i^Field = QE(N - Nd) (4) 

From this we can obtain the following expression for the electrophoretic mobility in absence 
of hydrodynamic interactions, where fi = 1/T is the mobility of a single monomer: 

£_ = N-N cl 
fio N + N cl 

For plotting © in figure @J we obtain Nqi by counting the average number of counter- ions 
found within 2<7o of the chain. The local force picture successfully describes the observed 
behaviour in absence of hydrodynamic interactions. 

In the linear response regime neither the chain structure nor the counterion distribution is 
affected by the presence of hydrodynamic interactions. Thus, the differences displayed in 
figures |3] and 0] can only be attributed to different frictional forces with the solvent acting on 
the PE-counterion compound. In the remainder of this paper, we will investigate this in more 
detail. 



3.4 Introducing hydrodynamic friction 

Figure 0] shows that the mobility strongly increases when hydrodynamic interactions are taken 
into account. In other words, the frictional forces the compound experiences from the solvent 
are reduced by the hydrodynamic interactions between the monomers and the counterions. 
The solvent interactions are no longer local, thus © and © have to be modified: 

^Solvent = -T cS (N,N C1 )v (6) 
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i 71 Field = Qeff(iV,iV CI ) J E (7) 

Here, r cff and Q cS are the a priori unknown effective friction and effective charge of the 
compound, that depend on the degree of polymerization N of the PE chain as well as on 
the number of co-moving counterions iV CI . Combining both equations leads to a general 
expression for the electrophoretic mobility: 

_ Q eS (N,N cl ) 
* T cS (N,N cl ) W 

Since the general functional dependence of r c g on N and Nci in the presence of hydro- 
dynamic interactions is non trivial, in all but limiting cases, © cannot be used to determine 
the electrophoretic mobility analytically. However, it can be applied to determine the effective 
friction from the electrophoretic mobility and the effective charge: 



eff 



Qeff/^. (9) 



3.5 Estimating the effective charge 

As before, the electrophoretic mobility is determined by Q, leaving the task to estimate the 
effective charge from the simulation. There are several possible ways to determine the number 
of counterions that shield the bare charge of the poly electrolyte. In section 13.31 we already 
used the straightforward way of counting the number of counterions within a tube of radius 
2a around the chain. 



Q { V = N-N cl (d<2a) (10) 

This method has the draw back, that the threshold, i.e. the size of the tube, is arbitrarily de- 
fined. Another frequently used way is the inflection criterion fl30u31[l . in which the integrated 
radial counterion distribution is calculated and the inflection point used as a threshold. 

Alternatively, we can rewrite (0) and use the electrophoretic mobility without hydrodynamic 
interactions to calculate the effective charge. 

Q$ = n(i-\^?) (ID 



i + A*r , 

This definition of the effective charge does not have any free parameters and since simula- 
tions without hydrodynamic interactions are computationally inexpensive, we can determine 
Q e s with high accuracy in this way. Further analysis showed that Q eS is insensitive to hydro- 
dynamic interactions in the linear response regime. 

In figure [5J we compare both estimators for the effective charge. For short chains, both 
estimators agree and coincide with the bare, unscreened charge of the poly electrolyte. In this 
regime, no counterion condensation is observed. 

For intermediate chains, the effective charge is reduced as it deviates from the bare charge 
and tends towards the Manning prediction. Here, the condensation parameter for the PSS 
system is £ = 3.12. In this regime, the simple estimator measures a higher effective 
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Figure 5: The effective charge Q cS as a function of chain length for both estimators and 
Q e g. The dotted line indicates the bare, unscreened charge of the polyelectrolyte, 
whereas the full line shows a prediction based on counterion condensation theory, 
with £ being the condensation parameter. 
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charge, i.e. not all condensed counterions that are included in the new estimator Q c{ { are taken 
into account. 

The so defined effective charge estimators and will now be used to quantify the 
effective friction of the polyelectrolyte-counterion compound. 




Figure 6: Effective friction of the poly electrolyte and the co-moving counterions as abtained 
from © using and . For comparison, the inverse diffusion coefficient is 
displayed, which is a measure for the hydrodynamic friction if the Einstein equation 
(fT2j) holds. The dashed lines indicate the approximate Debye screening length for 
the system. 



Figure [6] displays the effective friction of the compound in dependence of the length of the 
PE chain as obtained from © using the estimators and Q^j-. We compare the obtained 
result to the effective friction as it can be obtained via Einstein equation from measurements 
of the diffusion of the PE: 

r D = ^ (12) 
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Figure 7: The effective charge and effective friction per monomer as a function of the chain 
length shows an initially stronger decrease of the friction compared to the charge 
which leads to the observed increase in mobility. For long chains, both quantities 
approach a constant value (the lines are visual guides only). Again the dashed lines 
indicate the range of electrostatic screening. 
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For chains with contour lengths of the order of the Debye length, a clear deviation between 
r eff and T D is observed. Beyond this length scale, the effective friction of the PE increases 
almost linearly with chain length, where as the inverse diffusion scales with an exponent of 
m = 0.63 as shown in Section [3TT1 

In figure U\ we compare the effective friction and the effective charge per monomer. First, 
the effective friction shows a stronger decrease than the effective charge. This leads to the 
observed increase in mobility (see figure 0). For longer chains, both quantities approach a 
constant value. Again the relevant length scale seems to be the Debye length, beyond which 
both quantities become constant (i.e. linear scaling of Q c g and r ef j). 

This gives us a microscopic understanding of the effects, that lead to the experimentally 
observed length independent mobility of long flexible poly electrolytes. The hydrodynamic in- 
teractions between polyelectrolyte and counterions decrease the effective friction of the com- 
pound on length scales smaller than the electrostatic screening length. For compounds of 
longer chains, the hydrodynamic interactions are screened, and the friction increases linear 
with the size of the PE. 

As an application of the presented model, we will now look at the a possible way to increase 
the effective friction of the compound. 



3.7 Increasing the effective friction 



As suggested by Mayer et. al. 11321 13311 attaching a drag label ("molecular parachute") to one 
end of the PE, increases the hydrodynamic friction by a constant amount independent of the 
chain length. This leads to a changed scaling behaviour of the effective friction, which in turn 
can be seen in a restoration of the size dependent mobility even for long PE chains. 

In the simulations, the label consists of L = 40 uncharged monomers attached to one end 
of the PE chain. The hydrodynamic interactions with the fluid cause an additional force acting 
on the PE. Figure [8] demonstrates the effectiveness of this method. 

In the future, we will use the toolkit described in this article to analyse the ELFSE process 
in more detail. 



4 Conclusions 

We used a coarse-grained MD model to describe unlabeled and labeled free-solution elec- 
trophoresis of flexible poly electrolytes. The simulation results are in good quantitative agree- 
ment with experimental data. We used simulations without hydrodynamic interactions to ef- 
ficiently measure the effective charge of the polyelectrolyte as a function of the chain length. 
In a further step, we quantified the effective friction of the compound formed by the poylelec- 
trolyte and co-moving counterions. Our results show that on a length scale that is of the order 
of the Debye length for electrostatic screening the friction becomes linear in terms of the chain 
length. This increase in friction is exactly canceled out by the likewise linearly increasing ef- 
fective charge, leading to the well-known constant mobility for long flexible polyelectroyte 
chains. 
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Figure 8: The mobility of a long flexible polyelectolyte without label is independent of the 
chain length N (the line is a visual guide only). With a label of L = 40 uncharged 
beads attached to the polyelectrolyte, the size dependence of the electrophoretic mo- 
bility can be recovered. 
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